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ABSTRACT 

A  group  of  stations  in  the  North  American  Arctic  region  have  been 
analyzed  for  statistical  determination  of  temperatures  at  mandatory 
pressure  levels  p..   For  each  station  the  temperature  at  a  key  level 
(called  the  forcing-level  temperature)  peculiar  to  the  station  has 
been  forced  in  at  the  first  step,  and  retained  at  each  subsequent  step 
in  the  development  of  the  stepwise  regression  equations  giving  tem- 
perature at  the  mandatory  levels.   In  general,  eight-step  prediction 
equations  in  terms  of  other  temperatures  in  the  vertical  were  found  to 
give  specification  of  T(p  ),  with  percentage  explained  variance  of 
close  to  0.99,   As  a  result  of  this  definitive  property,  the  best- 
estimate  of  the  regional  atmosphere  which  is  conditionally  dependent 
upon  the  existence  of  an  extreme  1/^  probability  of  the  forcing-level 
temperature  is  obtained  with  a  high  degree  of  confidence. 
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1.   Introduction 

A  preliminary  set  of  model  atmospheres  giving  both  the  warm  and  cold 
extremes  at  certain  reference  elevations  with  a  1  ^  level  of  expectancy 
has  been  proposed  as  the  Preliminary  MIL-STD-210B  by  ETAC  (Environmental 
Technical  Applications  Center  of  the  USAF).   In  essence  this  model  of 
cold  and  warm  extreme  temperatures  is  based  upon  listing  the  coldest 
(warmest)temperatures  at  each  of  the  eight  standard  levels  listed  in 
Table  1.   These  extreme  temperatures  are  based  upon  observations  of  the 
location  of  extreme  temperature  noted  on  global  temperature  maps  at  the 
level  noted.   The  extremes  are  observed  to  occur  at  the  geographically 
diverse  locations  listed  in  Table  1. 

The  values  of  the  extremes  in  Table  1  were  computed  from  the  assumed 
existence  of  a  normal  distribution  of  temperature  in  both  the  coldest 
(January,  N.  Hemisphere)  and  warmest  (July,  N.  Hemisphere)  months.   Thus 
these  extremes  have  been  computed  based  upon  the  expectancy  (true  for  a 
normal  distribution)  of  finding  1^  of  all  temperature-departures  from  the 
sample  mean  at  each  location  and  level  at  2.3267  times  the  standard 
deviation  O   at  that  level  and  location.   In  the  process  of  making  the 
estimation  of  temperature-extremes,  no  allowance  has  been  made  for  the 
possibility  of  inter-level  correlations  which  have  been  shown  to  exist 
at  typical  stations  (Cole  and  Nee,  1965). 

Because  of  the  neglect  of  vertical  consistency  in  assembling  the  tem- 
perature-extreme data  for  the  Preliminary  MIL-STD-210B  climatic  model,  an 
unrealistic  range  of  temperatures  is  required  in  design  of  equipment 
needed  to  operate  in  these  atmospheres.   As  a  result,  it  has  been  proposed 
that  the  geographical  sites  which  are  proposed  in  Table  1,  as  correspond- 
ing to   world-wide  temperature  extremes  at  the  given  levels,  be  used  to 


TABLE  1.   Locations  of  proposed  extreme  temperatures  in  the 
cold  world-wide  and  warm  world-wide  cases  (after  ETAC). 


Co Id- Extreme 

Warm- Extreme 

Level 

Location 

January 
Mean 
(°C) 

STD. 
DEV. 
(°C) 

1^ 
Extreme 
(°C) 

Location 

July 
Mean 
(°C) 

STD. 
DEV. 
(°C) 

1^ 
Extreme 
(°C) 

SFC 

Oymyakon , 
USSR 

-50.6 

4.4 

-60.6 

Insalah, 
Algeria 

40.9 

3.9 

49.0 

850mb 

Oymyakon , 
USSR 

-35.1 

5.5 

-47.0 

Insalah, 
Algeria 

28.9 

2.2 

34.0 

CI 

700mb 

Hall  Beach, 

NWT 

-27.3 

6.5 

-42.4 

Babylon, 
Iraq 

17.0 

2.2 

22.1 

CI 

500mb 

Resolute, 

NWT 

-43.0 

4.6 

-53.2 

New  Delhi, 
India 

-4.3 

3.9 

4.0 

CI 

300mb 

Thule, 
Greenland 

-60.6 

2.2 

-66.0 

New  Delhi, 
India 

-25.8 

3.9 

-16.0 

CI 

200mb 

Thule, 
Greenland 

-59.7 

7.2 

-75.0 

Alert, 
NWT 

-42.3 

2.5 

-36.5 

150mb 

Karachi, 
Pakistan 

-65.0 

6.0 

-79.0 

Alert, 
NWT 

-43.3 

2.5 

-37.5 

lOOmb 

Singapore 

-83.0 

3.0 

-90.0 

- 

Thule, 
Greenland 

-43.8 

2.2 

-37.2 

Wl 
Wl 
Wl 


provide   "forcing  temperatures" at  the  level.   Based  upon  inter-level  tem- 
perature-statistics at  each  station,  the  most  likely  temperatures  at  the 
eight  additional  levels  (the  70  mb  level  is  to  be  included)  in  the  atmos- 
phere are  then  to  be  estimated  with  optimal  accuracy. 

It  is  to  be  noted  that  Table  2  lists  only  six  different  stations  in 
definition  of  the  cold-extreme,  and  five  in  definition  of  the  warm-extreme 
atmosphere.   Hence,  a  real-time  data  analysis  would  require  a  detailed 
climaticalogical  examination  of  11  different  stations  for  the  purpose  of 
establishing  regional  extremes, 
2.   A  statistical  model  for  regional  extreme  atmospheres. 

As  indicated  in  the  study  by  Cole  and  Nee  (1965),  high  values  of 
inter-level  simple  correlations  of  temperature  frequently  exist.   It  was 
decided  in  this  study  to  obtain  an  initial  data  file  on  magnetic  tape  for 
the  two-a-day  rawinsonde  stations  (all  in  the  North  American  sector) 
identified  by  code  sjTmbols  CI  and  Wl  in  Table  1,   These  symbols  are 
employed  to  suggest  the  cold  and  warm  atmosphere  feasibility  study.   The 
data  was  provided  by  the  Environmental  Data  Service,  National  Climatic 
Center,  through  the  financial  assistance  and  kind  cooperation  of  the  Com- 
mander, Naval  Weather  Service. 

The  other  stations  in  Table  1  not  bearing  the  code  symbols  CI  or  Wl 
have  also  been  provided  through  the  data  sources  just  listed.   However, 
these  secondary  data  records  are  not  available  in  the  convenient  summary 
form,  nor  in  such  abundant  population  samples  as  those  stations  which  were 
coded  CI  or  Wl  so  that  more  care  is  required  in  the  data-handling  of  these 
secondary  stations.   This  second  set  of  data-stations  is  the  subject  of 
ongoing  work  which  will  appear  as  a  later  report  in  this  series. 


For  the  stations  marked  code  symbols  CI  and  Wl,  a  stepwise  multiple 
regression  technique  was  employed.   This  program,  known  as  BIMED  02R,  is 
available  in  the  Program  Library  of  the  W.  R.  Church  Computer  Center  of 
the  Naval  Postgraduate  School,   The  temperature  data  was  arranged  in  a 
decreasing  pressure  sequence  but  with  standard  pressure-spacings  as  fol- 
lows (Table  2): 

TABLE  2. 


J=l 

J=2 

J=3 

J=4 

J=5 

J=6 

lOOOmb 

950 

9T)0 

850 

800 

750 

J  =  7 

700 

650 

600 

550 

500 

450 

J  =  12 

J  =  13 

400 

350 

300 

250 

200 

175 

J  =  18 

J  =  19 

150 

125 

100 

80 

70 

J  =  23 

Actually,  each  CI  and  Wl  radiosonde  was  listed  with  both  a  height  and 
temperature  at  each  of  these  levels.  Over  the  four  year  data  sample 
(1967-70),  these  stations  provided  the  following  population  samples: 

C-1  Stations  W-1  Stations 

Hall  Beach  (700mb)  243  samples  Alert  (200mb)  244  samples 

Resolute    (500mb)  237  samples  Alert  (150mb)  244  samples 

Thule  (300mb)  213  samples  Thule  (150mb)  236  samples 

Thule  (200mb)  213  samples 

The  program  BIMED  02R  permits  for  each  station,  the  specification  of 
a  different  dependent  variable  or  predictand  in  each  subproblem.   Also, 
the  number  of  independent  variables  to  be  tested  for  admission  to  regres- 
sion equation  can  be  limited.   In  the  problems  considered  here  an  upper 
limit  of  eight  independent  variables  was  found  to  be  given  optimal  speci- 
fication of  each  of  the  dependent  variables.   Another  convenient  feature 
of  the  BIMED  02R  program  is  that  the  forced  entry  into  the  regression  of 


the  temperature  value  being  tested  as  an  extreme  (the  so-called  "forcing 
temperature")  is  allowed,  even  though  its  contribution  to  the  explained 
variance  may  become  small.   For  example,  in  the  case  of  Resolute,  which 
corresponds  to  a  winter  extreme  at  500  mb,  the  multiple  regression  equa- 
tions at  the  various  "mandatory  levels"  were  set  up  in  the  forms: 

T^  =  A  +  B  T^^  +  C  X^  (1) 

Here  T^  is  the  8-level  predictand  chosen  alternately  as 

rri  ^      'Tp  ffi  rp  rn  rp  m  ryt  m 

^M  "  1*  ^4'  n*    nS*  ^17*   19'   21'  ^23 
where  the  subscript  M  is  the  temperature  identifier  associated  with  the 

pressure  level  indicated  in  Table  2.   In  summary,  the  predictands 

desired  are 

T, ,  the  temperature  at  1000  mb 

T, ,  the  temperature  at  850  mb 

T, ,  the  temperature  at  700  mb 

T^  ^ ,  the  temperature  at  500  mb 

T-c»  the  temperature  at  300  mb 

T. ,,  the  temperature  at  200  mb 

T^Q,  the  temperature  at  150  mb 

T-. ,  the  temperature  at  100  mb 

T  -,  the  temperature  at  70  mb 

The  form  X^  represents  here  a  linear  combination  of  up  to  7  independ- 
ent variables,  apart  from  the  forced  temperature  T.,  in  the  case  of 
Resolute,  NWT.   The  variable  X^  varies  from  case  to  case,  as  the  predictand 
is  changed  from  T^  to  T, ,  to  T-,  etc.  X^  usually  involves  temperature 
predictors  which  are  physically  close  to  the  level  being  specified  as 
the  dependent  variable  Z..   For  example,  at  Resolute,  T^  [=  T(850)] 
was  best-fit  by  the  following  expression: 


T,    =  .3745  -  .0849  T, ,  +  1.0  X, 
4  11       4 

X,  =  .3614T,  +  .85771^  -  .2892T,  +  .0819T, 
4        3        5        6        7 

+  .0902T„  +  .0598T,,-  .0486T,, 
9         14        17 

with  a  multiple  regression  coefficient  R  »  .9861.   In  general,  somewhat 

similar  results  applied  at  each  of  the  other  mandatory  pressure  levels 

1000,  700,  500,  300,  200,  150,  100,  70  mb  except  that  the  multiple  R 

expressing  T  [=  T(IOOO)]  and  T.^[=  T(300)]  in  terms  of  T^,  and  X^  were 

somewhat  lower  than  the  remaining  correlations  in  the  set.   Similar 

statements  applied  to  all  stations  at  both  extreme- times  January  and 

July, 

The  function  X^  is  called  in  this  study  a  "structure- function"  in 

that  to  a  large  degree  it  accounts  for  most  of  the  variation  of  the  T., 

M 

being  predicted,  even  though  the  part  played  by  the  forcing  variable 
is  by  no  means  negligible. 

A  least  squares  fit  of  form  Eq.  (1)  may  be  shown  to  produce  residuals 
(computed  against  the  observed  values)  which  are  distributed  normally 
relative  to  the  regression-value,  together  with  a  variance  given  by  the 
standard  error  of  estimate  a  =  a  Jn.  -  R^  .   In  the  set  of  dependent 
variable  regressions  tested  here,  nxjst  of  the  multiple  regression  coeffi- 
cients had  values  well  in  excess  of  .99  (with  the  exception  of  the  two 
levels  already  noted,  1000  mb  and  300  mb) ,  so  that  a  near  functional 
relationship  existed  between  the  predictor  and  the  predictand  set. 
The  specification  statistics  for  the  different  stations  in  the  North 
American  sector  of  Table  1  are  listed  in  Section  3. 


3.   Regression-specification  of  mandatory- level  temperatures  using  the 
forcing-level  temperature. 

The  stepwise  regression  procedure  of  BIMED  02R  was  employed  with  the 
appropriate  forcing-level  temperature  required  at  the  first  step  of  the 
regression  and  retained  in  the  selection  of  the  following  seven  predictors. 
These  seven  predictors  were  allowed  to  enter  the  stepwise  regression  in 
accordance  with  the  programmed  requirement  that  at  each  step  k,  the  variable 
added  explained  a  maximum  of  the  unexplained  variance  remaining  after  (k-1) 
selections  (k=2, ....... ,8) . 

Tables  3a,  4a,..., 9a  show  the  results  of  the  stepwise  regression  applied 
to  the  specification  of  the  temperatures  T,,  T,,  T,,  T  , ,  T,^*  T^,,  T,q, 
Tp, ,  T-»  at  mandatory  levels  by  an  equation  of  form  (1),  using  an  eight- 
predictor  equation  for  each  T..  just  listed.   Table  3a  lists  for  the  January 
full  data-sample  at  Hall  Beach  the  mean,  standard  deviation,  multiple  cor- 
relation coefficient  and  standard  error  a^at   each  level  other  than  at 
forcing  level.   Table  3b  lists  analogous  results  for  a  nominal  10^  of  the 
January  data  sample  corresponding  to  the  10^  cold  extreme  sample. 

Table  4a  and  4b  list  similar  results  for  Resolute  in  January  with 
T, ^[=  T(500)]  as  the  forced  variable.   Tables  5  and  6  list  analogous 

results  for  Thule,  using  T   and  T^  as  forcing  variables,  respectively. 

2 
It  is  to  be  noted  that  the  fractional  explained  variance  R  of  the  predic- 

tand  in  each  case,  does  not  fall  of  appreciably,  in  general,  as  one  pro- 
ceeds along  any  row  of  these  tables  from  part  (a)  to  part  (b) . 

Tables  7,  8,  9  similarly  list  the  July-warm  statistics,  using  respec- 
tively the  full  data-samples  at  Alert  [with  forcing  temperature  T  =  T(200)], 
Alert  [T  =  T(150)]  and  Thule  [T  =  T(IOO)]  in  the  left  half  of  the  tables. 

J  «J 

The  prediction  results  corresponding  to  the  nominal  10^  of  the  warmest 
July  forcing-level  temperatures  appear  on  the  right  half  of  these  tables. 


The  selection  of  nominal  10%  extreme  cases  was  made  by  separating  out  of 
the  full  sample  all  cases  where  the  forcing  temperature  T  lay  in  the  ranges 
T  <  T  -  1.2817a  for  cold  extremes 
Tj  >  Tj.  +  1.2817a  for  warm  extremes 
The  criteria  used  in  establishing  these  90/^  nominal  extremes  was  that  the 
full  data  samples  were  normally  distributed,  a  valid  first  approximation. 
The  actual  number  of  extreme  data  cases  is  listed  at  the  top  of  Tables  3b, 
....,  9b,  for  each  case  under  consideration. 

In  order  to  determine  the  multiple  correlation  coefficients  and  standard 
errors  for  the  right-side  10^  extreme  samples  of  Tables  3b,..., 9b,  a  regres- 
sion equation  of  identical  form  to  that  developed  for  the  full  data-base 
predictands  was  generated  for  testing  the  lOfo   extreme-case  data  samples. 
The  test-predictors  for  the  extreme-data  were  then  considered  to  be  the 
forcing  variable  T  and  the  structure-variable  X  of  form  determined  for 
the  full-sample  cases.   The  high  values  of  the  multiple  correlations  thus 
determined  for  the  January  extreme-data  cases  serves  to  verify  the  concept 
that  the  predictands  corresponding  to  the  extreme-data  may  be  anticipated 
through  a  joint  knowledge  of  the  forcing  variable  T  and  of  the  structure 
function.   The  somewhat  smaller  correlation  coefficients  for  the  July  data 
is  due  to  the  much  smaller  values  of  variance  in  summer,  particularly  in 
the  lower  stratosphere. 

A  listing  of  the  nominal  extremes  of  the  forcing  temperature  for  other 
probability  levels  of  interest,  e.g.,  the  5^,  the  1^  was  also  made  assuming 
the  existence  of  a  Gaussian  distribution  of  temperatura  and  appears  in 
Table  10.   A  separation  of  the  10^  nominal  extremes  of  T  was  then  made 
wherever  possible  into  the  classes  suggested  in  Table  10,  that  is  beyond 
the  cOl  probability  level;  in  the  probability  range  .05  to  .01  and  in  the 


10 


range  .10  to  .05.   The  results  of  these  probability  stratifications  on  the 
predictands  at  T- ,  T,,  T^,  T   ,  T^>  T^,  T   ,  T   are  discussed  in  Section  4, 
These  additional  critical  extremes  are  given  by 

T  (.05)  =  T  +  1.640  a  (2) 

*J  \J  «J 

T  (.01)  =  T  +  2.3267a  (3) 

J  J  J 

Equations  2  and  3  tend  to  give  excessive  estiimtes  of  the  critical 
extremes  when  the  observed  O     is  too  large  to  be  representative  of  a  normal 
distribution.   This  happens  primarily  in  the  winter  season  at  low  strato- 
spheric levels,  e.g.,  Thule  at  300  mb.   At  this  locale  there  tends  to  be 
bimodal  alternations  between  warm  and  cold  stratospheric  regimes  and  an 
unrepresentatively  large  o     results  from  this  sample.   The  ETAC  estimate 
(Table  1)  of  the  lf>   extreme  temperature  at  Thule  at  300  mb  in  winter  also 
seems  to  reflect  this  same  difficulty. 

The  following  temperature  classes  were  then  considered  as  input  data  at 
the  CI  and  Wl  station  forcing  levels,  respectively 
Cold  classes:   (a)  T^<   T^^^;  (b)  T  q^<  T^<  T  ^3;  (c)  T  ^3  <  T^<  T^^^ 

Warm  classes:   (a)  ^>T^q^;  (b)  \o-^>  Tj>  T ^^^;    (c)  \o^  >  ^j>  \io 

In  most  cases  studied,  reasonable  sample-sizes  existed  in  each  class  (a), 
(b),  (c).   However,  in  several  cases,  the  extremes  did  not  reach  the  1^ 
nominal-value  predicted  from  an  assumed  Gaussian  distribution. 
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TABLE  10.   Nominal  temperature-extremes  at  the  indicated 
stations  and  pressure  levels  at  the  10^,  5^  and  1^  extreme 
values  according  to  a  Gaussian  distribution. 


90^  extreme 

Station  T  ^_  T  _  T  ^^  (°C)   Forcing  Level  T    sample  size 

Hall  Beach  -35.6  -37.9  -42.3 

Resolute  -47.1  -48.5  -51.3 

Thule,  Wi  -62.0  -63.8  -67.3 

Thule,  Wi  -63.5  -66.0  -70.6 

Thule,  Su  -41.3  -40.8  -39.8 

Alert,  Su  -39.7  -39.1  -37.9 

Alert,  Su  -38.6  -37.6  -35.7 


T     =  T(700) 

vJ 

26 

T     =  T(500) 

25 

T     =  T(300) 

13 

T     =  T(200) 

J 

18 

T      =  T(IOO) 

18 

T     =  T(150) 

22 

T     =  T(200) 

J 

18 

19 


4.   Temperature-estimates  at  mandatory  levels  corresponding  to  1/^  extremes 
of  forcing— level  temperatures. 

a.   Winter  extreme  atmospheres 

Corresponding  to  the  four  winter  extreme  cases,  the  multiple  regres- 
sion methods  which  led  to  Tables  3,  4,  5,  6  also  lead  to  the  results  of 
Tables  11,  12,  13,  14.   In  the  latter  tables,  the  nominal  10^  extreme  sample 
is  decomposed  into  the  subsample  extreme-classes  (a),  (b) ,  (c).   These 
classes  were  defined  in  Section  3,  but  their  definition  is  also  implicitly 
given  in  the  top  row  of  each  table  (e.g.,  Table  11:  (a)  T  q- ;  (b)  T  _.  - 
T  ^j.;  (c)  T  ^c  ~  T  ^_  .   Corresponding  to  the  stratified  data-sample  at  each 
forcing  level  there  exist  also  conditional  data-sets  at  each  of  the  manda- 
tory levels  above  and  below  the  forcing  level.   In  virtually  all  cases, 
the  mandatory- level  data  sets  in  classes  (a),  (b),  (c)  exhibit  well-defined 
trends  between  the  class-mean  temperatures  at  all  levels.   This  discovery 
afforded  credibility  to  the  listed  extreme  values  of  the  class  (a)  mandatory- 
level  results  of  Tables  11,  12,  14.   Note  that  in  each  stratification  (a), 
(b),  (c) ,  the  same  sample  size  exists  at  each  level  within  each  class  re- 
gardless of  the  level  under  consideration. 

As  noted  in  Section  3,  no  nominal  1^  data  class  was  realized  in 

Table  13;  however,  the  class  (b)  sample  T  =  T(300)  comprised  a  de  facto 

extreme  1^  data-sample  by  actual  count.   In  the  comparison  of  T,  with  T  at 

b       c 

mandatory  levels  over  Thule  (Table  13),     there    existed  consistent 
Inter-class  temperature  differences,  again  at  all  levels.   In  all  cases, 
the  reason  for  these  consistent  inter-class  temperature  deviations  is  that 
the  mandatory  level  temperatures  have  been  found  to  be  strongly  correlated 
with  the  forcing -levels  temperature  through  the  multiple  regression 
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Level 

Teicperature-means    (   C) 

Std. 

deviations 

(°C) 

(mb) 

in  classes 

relative   to   class 

-mean 

Ca) 

(b) 

(c) 

^.01 

T        -T 
.01    \05 

T        -T 
\05    \lO 

a 
a 

^b 

a 
c 

1000 

-35.24 

-35.22 

-34.62 

2.588 

4.047 

3.876 

850 

-39.52 

-35.71 

-32.34 

1.421 

2.475 

i.943 

700 

-44.96 

-39.55 

-36.60 

1.746 

1.175 

0.799 

500 

-44.04 

-45.90 

-47.32 

0.948 

3.363 

1.233 

300 

-41.02 

-45.95 

-47.41 

0.937 

5.472 

3.315 

200 

-42.96 

-46.22 

-47.16 

0.924 

6.064 

4.327 

150 

-45.20 

-47.32 

-48.91 

0.767 

7.804 

6.220 

100 

-48.64 

-51.40 

-52.25 

1.284 

8.999 

6.849 

70 

-50.28 

-55.47 

-55.42 

2.047 

10.185 

6.527 

Number 

in 

class 

5 

13 

8 

TABLE  11.   Temperature  means  at  mandatory  levels  at  Hall  Beach  (winter), 
corresponding  to  temperature  extremes  of  Tj(700")  at  the  (a)  I'p   class, 
(b)    1-5^  class  and  (c)  5-IO5U  class  of  probability.   The  standard  devia- 
tions 0,0,0     of  the  observed  temperatures  within  each  class  are  shown 
in  the  ^   ^   ^  right  half  of  the  table. 
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■40        -20 
Temp.rO 


FIG.  1.   The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Hall  Beach,  NWT,  corresponding 
to  the  class  (a)  set  of  l'^   cold  extreme  occurrences  of  T  (700)  of 
Table  11.   The  thin  solid  line  depicts  the  January  mean  Hall  Beach 
vertical  sounding  (1967-70).   The  dashed  curves  the  world-wide  Pre- 
liminary MIL-STD  210  extreme  atmospheres. 
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Level 

Temperature-means    (   C) 

Std. 

deviations 

(°C) 

(mb) 

in  classes 

relative  to  class 

-mean 

(a) 

(b) 

(c) 

^.01 

T         -T 
\01   \05 

T        -T 
.05      .10 

a 

a 

^b 

O 
c 

1000 

-39.00 

-38.77 

-33.31 

0.748 

2.125 

3.718 

850 

-35.50 

-31.57 

-30.14 

0.698 

1.463 

3.038 

700 

-40.37 

-36.43 

-34.45 

0.742 

1.308 

2.274 

500* 

-51.73 

-49.70 

-47.59 

0.173 

0.782 

0.342 

300 

-54.03 

-56.20 

-60.48 

1.225 

2.027 

4.070 

200 

-55.20 

-55.24 

-59.15 

0.648 

2.280 

3.561 

150 

-59.37 

-58.09 

-61.12 

0.881 

2.489 

3.098 

100 

-65.70 

-61.00 

-64.59 

0.860 

3.745 

3.248 

70 

-70.77 

-62.50 

-65.88 

0.480 

5.748 

4.064 

Number 

in 

class 

3 

7 

15 

TABLE  12.   Temperature  means  at  mandatory  levels  at  Resolute  (winter), 
corresponding  to  temperature  extremes  of  T  (500)  at  the  (a)  1%   class, 
(b)  1-5^  class  and  (c)    5-10^1  class  of  probability.   The  standard  devia- 
tions C7  ,  a.  ,  a  of  the  observed  temperatures  within  each  class  are 
shown  ^       ^  in  the  right  half  of  the  table. 
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-20 
Temp.(  C) 

FIG.  2.   The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Resolute,  NVTT,  corresponding  to 
the  class  (a)  set  of  1^  cold  extreme  occurrences  of  T^(500)  of 
Table  12.   The  thin  solid  line  depicts  the  January  mean  Resolute 
sounding  (1967-70).  The  dashed  curves  are  Identical  to  those  of 
Figure  1. 
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1000 


relationships  summarized  in  Tables  3  through  9.  These  relationships  also 
apply  across  classes  (a),  (b),  (c),  and  give  rise  to  the  inter-class  dif- 
ferences found  in  Tables  11,...,  17. 

The  listings  in  Tables  11,  12,  13,  14  of  the  class  (a)  set  of 
temperatures  [or  the  (b) -class  if  class  (a)  is  non-existen^  at  mandatory 
levels  makes  possible  the  presentation  of  the  results  in  graphical  form. 
Thus  the  expected  (cold)  temperature  data  corresponding  to  the  1/S  forcing 
temperature  appropriate  to  Tables  11,  12,  13  and  14  have  been  reproduced 
as  Figs.  1,  2,  3,  4.   Here  the  heavy  solid  line  depicts  the  cold  or  ly^ 
"extreme"atino8phere  in  each  case,  whereas  the  thin  solid  line  depicts  the 
January  (1967-70)  mean  sounding  in  each  case.   Note  that  the  temperature 
deviations  between  "mean"  and  "extreme"  in  Figs.  1,  2,  3,  4  are  indicative  also 
of  the  class  (a)  to  (c)  deviations  at  the  mandatory  levels  of  Tables  11, 
12,  13,  14. 

b.   Warm  extreme  atmospheric  cases. 

The  procedure  in  these  cases  consists  of  an  analysis  of  the  nominal 
10^  extreme-warm  July  cases  at  Alert  (  [Tj  =  T(200)],  [T^  =  T(150)]j   and  at 
Thule  [T  =  T(IOO)]  after  these  extremes  have  been  subdivided  into  the 
nominal  probability  classes  (a),  (b),  (c),  similar  to  the  procedure  of 
Section  4(a). 

In  two  of  the  three  cases,  the  class  (a)  nominal  probability  T^< 
T  ^.  did  not  occur  (see  Tables  15  and  17),  but  an  actual  case  count  reveals 
that  in  either  case  there  were  2  to  4  data-samples  in  class  (b).   These 
samples  were  taken  as  representative  of  the  actual  1^  extreme  warm  situa- 
tion at  the  respective  forcing  levels,  and  also  at  the  regression- 
dependent  mandatory  levels. 
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Level 

Temperature-means    (   C) 

Std.    deviations 

(°C) 

(mb) 

in  classes 

relative  to  class 

!-mean 

(b) 

(c) 

T        -T 
^.01    \05 

T         -T 
\05   \lO 

% 

O 
c 

k— — — — — 

1000 

-25.20 

-20.27 

2.371 

4.942 

850 

-24.83 

-21.87 

1.637 

3.736 

700 

-29.73 

-28.36 

5.287 

3.858 

500 

-42.90 

-42.26 

2.032 

3.822 

300 

-64.47 

-62.60 

.173 

.341 

200 

-62.83 

-61.60 

.526 

2.294 

150 

-63.73 

-62.80 

.911 

2.889 

100 

-67.93 

-67.38 

1.496 

2.581 

70 

-71.50 

-70.29 

2.921 

3.026 

Number 

in 

class 

3 

10 

TABLE  13.   Temperature  means  at  mandatory  levels  at  Thule,  Greenland 
(winter),  corresponding  to  temperature  extremes  of  T  (300)  at  the  (b)  1-5^ 
class,  and  (c)  5-10^^  class  of  probability.   The  population  of  class  (a)  was 
zero.   The  standard  deviations  o     and  o     of  observed  temperatures  in  the  (b) 
and  (c)  classes  are  shown  in  the  right  side  of  the  table. 
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FIG.  3.   The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Thule,  Greenland,  corresponding 
to  the  1^  cold  extreme  set  of  forcing  temperatures  T  (300)  of  Table 
13.   The  thin  solid  line  depicts  the  January  mean  Thule  sounding 
(1967-70).  The  dashed  curves  are  identical  to  those  of  Figure  1. 
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Level 

(mb) 

o 
Tecr.per.Tture-nieans    (    C) 

in   classes 

Std.    deviations    (°C) 
relative   to   class-mean 

(a) 

T 
\01 

(b) 

T         -T 
\01      .05 

(c) 

T        -T 
\05       .10 

a 
a 

^b 

a 
c 

1000 
850 
700 
500 

300 

* 

200 

150 

100 

70 

-9.87 
-5.30 
-12.67 
-29.97 
-57.03 
-74.70 
-73.57 
-73.43 
-77.50 

-12.17 
-15.31 
-19.73 
-34.99 
-59.19 
-68.61 
-66.81 
-71.30 
-75.46 

-14.51 
-17.09 
-23.68 
-38.21 
-60.40 
-64.80 
-63.74 
-68.29 
-71.52 

4.302 
2.269 
0.806 
1.304 
0.695 
.084 
0.173 
0.733 
1.145 

3.849 
3.179 
1.608 
2.732 
1.041 

1.607 
2.496 
2.413 
3.156 

8.440 
5.024 
4.870 
3.773 
2.115 
.480 
2.216 
2.339 
2.767 

Number 

In 
class 

3 

7 

8 

TABLE  14.   Temperature  means  at  mandatory  levels  at  Thule  (winter),  corres- 
ponding to  temperature  extremes  of  Tj(200)  at  the  (a)  1^  class,  (b)  1-5^ 
class  and  (c)  5-10^^  class  of  probability.   The  standard  deviations  O   ,  a  , 
O  of  the  observed  temperatures  within  each  class  are  shown  in  the  right 
half  of  the  table. 
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Temp.(  C) 
FIG.  4.   The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Thule,  Greenland,  corresponding 
to  the  15^  cold  extreme  occurrences  of  T  (200)  of  Table  14.  The  thin 
solid  line  depicts  the  January  mean  Thule  sounding  (1967-70) .  The 
dashed  lines  are  identical  to  those  of  Figure  1. 
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Level 

Temperature-means  (  C) 

Std.  deviations 

(°C) 

(mb) 

in  classes 

relative  to  class 

-mean 

(b) 

(c) 

-^.01-^.05 

T    -T 
\05  \lO 

^b 

a 
c 

1000 

0.05 

1.48 

.350 

1.147 

850 

-6.65 

-4.07 

.650 

2.086 

700 

-16.15 

-13.09 

.250 

1.953 

500 

-31.30 

-27.79 

2.000 

4.086 

300 

-40.05 

-45.54 

2.250 

3.443 

200* 

-36.70 

-38.31 

.600 

.330 

150 

-38.90 

-39.67 

1.000 

.575 

100 

-40.95 

-40.47 

.350 

.860 

70 

-40.25 

-39.56 

1.150 

1.498 

Number 

in 

class 

2 

16 

TABLE  15.   Temperature  means  at  mandatory  levels  at  Alert  (summer),  corres- 
ponding to  temperature  extremes  of  T1200)   at  the  (b)  1-5/^  class  and  (c)  5- 
105^  class  of  probability.   The  population  of  class  (a)  was  zero.   The 
standard  deviations  Q.  >  ^  of  the  observed  temperatures  within  classes  b^ 
and  £  are  shown  on  the  right  side  of  the  table. 
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FIG.  5.  The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Alert,  NWT,  corresponding  to  the 
1^  warm  extreme  occurrences  of  T  (200)  of  Table  15.  The  thin  solid 
line  depicts  the  July  mean  Alert  sounding  (1967-70).   The  dashed 
curves  are  identical  to  those  of  Figure  1. 
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Level 

Temperature-means    (    C) 

Srd. 

deviations 

(°C) 

(mb) 

In  classes 

relatl 

ve   to   class 

-mean 

(a) 

(b) 

(c) 

^.01 

T         -T 
^01      .05 

T        -T 
\05      .10 

0 
a 

^b 

a 
c 

1000 

0.00 

2.68 

2.29 

.030 

1.384 

2.734 

850 

-4.65 

-3.93 

-3.21 

2.650 

2.318 

2.067 

700 

-13.75 

-12.48 

-11.94 

2.650 

2.467 

2.689 

500 

-26.25 

-26.42 

-26.48 

3.050 

1.117 

2.317 

300 

-43.50 

-49.23 

-46.78 

1.200 

1.168 

3.177 

200 

-38.10 

-38.83 

-39.09 

.800 

.619 

.909 

150* 

-38.20 

-38.80 

-39.57 

.300 

.224 

.226 

100 

-40.60 

-40.33 

-40.33 

.000 

.500 

.469 

70 

-39.45 

-39.93 

-39.43 

.349 

.472 

1.471 

Number 

in 

class 

2 

6 

14 

TABLE  16.   Temperature  means  at  mandatory  levels  at  Alert  (summer),  corres- 
ponding to  temperature  extremes  of  T/150)  at  the  (a)  l/o  class,  (b)  1-5^ 
class  and  (c)  5-10^^  class  of  probability.   The  standard  deviations  O  ,    O  , 
O     of  the  observed  temperatures  within  each  class  are  shown  in  the 
right  side  of  the  table. 
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FIG.  6.  The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Alert,  NWT,  corresponding  to  the 
1^  warm  extreme  occurrences  of  T-(150)  of  Table  16 o  The  thin  solid 
line  depicts  the  July  mean  Alert  sounding  (1967-70).  The  dashed 
curves  are  identical  to  those  of  Figure  1. 
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Level 

Temperature-means  (  C) 

Std.  deviations 

(°C) 

(mh) 

in  classes 

relative  to  class 

-mean 

(b) 

T    -T 
\01   .05 

(c) 

T    -T 
.05   .10 

% 

a 

c 

1000 

0.90 

2.56 

.100 

2.821 

850 

-2.42 

-3.72 

3.544 

1.943 

700 

-10.10 

-11.43 

3.412 

3.283 

500 

-25.80 

-26.69 

4.513 

2.884 

300 

-45.25 

-46.81 

2.296 

2.772 

200 

-39.80 

-39.85 

1.404 

1.211 

150 

-39.80 

-40.86 

1.931 

.875 

* 

100 

-40.20 

-41.16 

.122 

.105 

70 

-39.92 

-40.82 

1.803 

.699 

Number 

in 
class 

4 

14 

TABLE  17,  Temperature  means  at  mandatory  levels  at  Thule  (summer),  corres- 
ponding to  temperature  extremes  of  Tj(100)  at  the  (b)  1-5^  class  and  (c) 
5-10>  class  of  probability.   The  population  of  class  (a)  was  zero.   The 
standard  deviations  O^y    O     of  the  observed  temperatures  within  classes  b^ 
and  £  are  shown  on  the  right  side  of  the  table. 
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FIG.  7.   The  heavy  solid  line  shows  the  mean  regression-determined 
vertical  temperature  sounding  over  Thule,  Greenland,  corresponding 
to  the  1^  warm  extreme  occurrences  of  Tj(lOO)  of  Table  17.  The  thin 
solid  line  depicts  the  July  mean  Thule  sounding  (1967-70).   The  dashed 
curves  are  identical  to  those  of  Figure  1. 
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Table  15  thow*  that  class  (b)  waxm  atmosphere  Is  warmer  than  the 
class  (c)  case  in  the  layer  AOO  mb  to  123  mb.   Above  125  mb  Fig.  5  shows 
that  the  deviation  between  the  mean  case  and  the   listed  extreme  at  any 
level  Is  small.   For  example,  the  standard  deviation  a  at  Alert  at  the 
100  mb  level  Is  1.26°C  (and  is  smaller  still  at  70  mb).   Below  400  mb,  the 
temperature  of  the  class  (b)  extreme  atmosphere  is  consistently  colder 
than  that  of  the  class  (b)  atmosphere.    All  of  these  features 
are  also  to  be  observed  in  Fig.  15,  which  gives  a  comparison  of  the  1/^ 
warm  extreme  ataosphere  and  the  July  mean  (1967-70)  over  Alert. 

A  similar  set  of  conclusions  applies  to  the  regression-dependent 
class  (a)  ataospharas  over  Alert  based  upon  T  =  T(150)  as  forcing  variable, 
Here  the  results  included  in  Table  16  indicate  that  the  500  mb  tempera- 
ture is  a  level  of  saall  interdlumal  variability,  and  is  in  fact  the 
location  of  a  crossover  between  the  T(extreme)  and  T(mean)  curves.   Below 
500  mb,  the  inferred  temperature  structure  corresponding  to  a  forcing 
temperature  at  the  1^  warm  extreme  at  p  =  150  mb  shows  an  atmosphere 
slightly  colder  than  the  July  mean. 

The  sujuaer  extreme  atmosphere  at  Thule,  based  upon  a  forcing- level 
warm  extreme  In  T.  ■  T(IOO)  is  listed  in  Table  17  and  also  graphically 
depicted  in  Fig.  7.   The  results  summarized  in  column  (b)  of  Table  17 
show  that  the  extrcae  warm  case  at  100  mb  is  associated  with  mandatory 
levels  which  are  consistently  warmer  at  all  levels  (other  than  the  surface) 
than  in  the  class  (c)  atmosphere. 
5 .   Recommendat  ions 

For  the  best  estiaate  of  the  forced-level  extreme  atmospheres  con- 
sidered here,  one  merely  reads  off  the  temperature  of  the  left-hand  extreme 
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class  of  Tables  11,..., 17.   To  be  more  realistic,  one  should  allow  a  tem- 
perature-range of  one  standard-deviation  ±  0     about  any  class  (a)  extreme 

cl 

atmosphere.   However,  if  a  single  temperature-estimate  must  be  used  at  each 

level,  the  mean  T  -profile  as  a  function  of  P  is  to  be  recommended.   If 
a 

class  (a)  is  lacking,  then  the  T, -profile  should  be  used. 

The  work  presented  here  does  not  include  certain  other  levels  listed 
in  the  preliminary  MIL  STD  210-B  atmosphere.  The  four  years  of  climatologi- 
cal  data  for  the  years  1967-70  for  these  additional  stations  (see  Table  2) 
is  still  undergoing  data  analysis.  The  results  for  OJMJAKON,  Siberia,  have 
been  analyzed  statistically  and  have  been  found  to  determine  regression- 
generated  mandatory- level  temperatures  much  like  those  derived  for  the 
North  American  stations.   It  is  expected  that  a  study  of  the  additional 
stations  listed  in  Table  2  will  yield  fruitful  results. 
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APPENDIX 

It  is  not  feasible  due  to  the  limited  sample  sizes  to  employ  a  t- 
test  on  the  difference  of  means  for  the  1/^  class  relative  to  the  1-5^ 
class  of  Tables  11,..., 17,  However,  a  well-defined  test  of  signifi- 
cance of  class-means  is  afforded  by  the  t-test  applied  to  the  differ- 
ence of  sample-means  obtained  by  the  results  of  the  predictions  of 
the  neans  of  the  10/^  extreme  relative  to  the  full  sample  predictions 
(e.g.,  the  means  contained  in  Tables  3,  4,.. .,9). 

The  t-test  on  the  predictions  of  the  full-sample  and  the  10^  ex- 
treme involves  at  every  mandatory  level  and  station  a  test  for  signi- 
ficance of  the  pooled  t-statlstic  after  prediction. 

Here  the  T's  are  the  sample  temperature  means  at  the  mandatory  levels 

based  upon  the  identical  prediction  equations.  These  means  are  listed 

in  columns  2  and  7  of  Tables  3,  4,..., 9.  Q  is  the  standard  error  of 

iu 

estimate  listed  in  column  5,   and  0_     ,n  ^s  that  of  colcfiui  10,   associated 

£,  .  lU 

with  the  prediction  of  the  extreme  sample.  The  number  of  degrees  of 

freedom  associated  with  these  o  -values  are  N-9  and  n-2  respectively, 

where 

N  =  full-sample  size 

n  =  nominal  10^  extreme  sample  size 

The  pooled  tp-parameter  is  then  tested  using  frjf  +  n-ly  degrees  of  freedom 

for  tp. 
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At  the  70-nib  level  in  Table  3,  the  t  -statistic  based  upon  the 
listed  difference  of  the  means  is  t  =  3.113.   Using  the  number  of 
degrees  of  freedom  N  +  n-11  =  258,  the  sample  t  -statistic  is  signi- 
ficant at  a  level  of  99. 9^^  probability. 

In  a  similar  manner  all  sets  of  differences  of  full-sample  and 
"extreme"  sample  means  obtained  by  the  prediction  method  of  this 
study  prove  to  be  significant  at  levels  in  excess  of  99/^  confidence. 
This  may  be  verified  by  applying  all  the  statistics  of  Tables  3,  4, 
...,9(a,b)  to  the  test-statistic  t  of  Eq,  (A-1). 

In  this  Appendix,  t-test  procedures  for  considering  definitive 
tests  of  significant  differences  of  means  between  full-  and  extreme- 
sample  statistics  at  mandatory  levels  have  been  set  forth.   In 
defining  the  extreme-sample,  the  10^  probability  extreme  at  the 
forcing  level  has  been  considered  as  the  basis  of  the  sampl^  includ- 
ing the  resulting  regression-statistics.  Beyond  the  10^  probability 
level,  reduced  sample  sizes  preclude  the  usual  t-test  procedures. 
However,  the  stratification5 based  upon  extreme  sub-samples  (a),  (b), 
(c)  of  Tables  11,  12,..., 17  give  rise  to  regression-generated 

T  (meani  which  might  well  have  been  extrapolated  from  the  10^  sub- 
a 

sample  results  of  Tables  3,  4,..., 9,  respectively. 
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